function y = e_qp1(x)
% Extended Quadratic Penalty 1 Function   
%x0=[1,1,...1].  
%fmin=n/2
%Author:shi yangyang
n=length(x);
u=ones(1,n-1);
v=ones(1,n);
for i=1:1:n-1
    u(i)=(x(i)*x(i)-2)^2;
end
for i=1:1:n
    v(i)=x(i)*x(i);
end
y=sum(u)+(sum(v)-0.5)^2;